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1. INTRODUCTION 

It is generally expected that two-flavor QCD 
undergoes a high temperature chiral symmetry 
restoring phase transition at zero quark mass, 
with 0(4) critical behavior [jjj. Verifying this ex- 
pectation is important for understanding the phe- 
nomenology of the transition and for facilitating 
an extrapolation of simulation data to physical 
quark masses. Lattice QCD at fixed Nt with stag- 
gered fermions is similarly expected to exhibit at 
least 0(2) universality, with 0(4) emerging in the 
continuum limit. For present purposes we do not 
distinguish 0(2), 0(4), or mean field critical be- 
havior. 

One test of universality compares critical ex- 
ponents. Comparing the critical scaling function 
gives further insight. We use the standard corre- 
spondence between QCD variables and O(N) spin 
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variables, which identifies quark mass m q /T with 
magnetic field h, inverse gauge coupling 6/g 2 with 
temperature T/T c (0), and chiral condensate •tpt/j 
with magnetization M. A critical point is ex- 
pected to occur at zero quark mass and nonzero 
coupling 6/g 2 (0). For studies at fixed Nt we use[0 

h = am q N t 

t = T/T c (0) - 1 - 6/g 2 - 6/g 2 (0) 

Critical scaling theory predicts that for small 
quark masses we have the scaling relation (|| 

■ipiph~ 1/s = /qcd{x) = c y fo(4)(c x x) 

where x — th~ x ^ s , c x and c y are scale constants, 
Jqcd(x) is the critical scaling function for QCD 
and fo(4)( x ) is that for 0(4). Not shown are non- 
leading nonscaling contributions analytic in t and 
h. 

2. A CRITICAL EXPONENT 

We extended an old data sample at Nt = 4, 
which was done with the standard one-plaquette 
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gauge action plus two-flavor staggered fermion ac- 
tion. Our new sample decreases the quark mass 
to am q — 0.008 and increases the spatial lat- 
tice size L to 24. We also reanalyze old data at 
N t = 6, 8, 12 f§. The latter all have L = 2N t . 

The mixed plaquette/chiral condensate suscep- 
tibility 

Xmt = d(^)/d(6/g 2 ) 

has a peak, which can be used to locate the 
crossover. We extrapolate the peak height to in- 
finite volume as shown below. 
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The infinite volume value should scale as 

Xm t) max ~ (am,)^ 1 )/^ « (am,)- ' 33 

The figure below shows a strong departure from 
scaling at present quark masses, in agreement 
with similar findings of the JLQCD and Bielefeld 
groups for a variety of susceptibilities ||]. 
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3. SCALING FUNCTION AT N t = 4 

To find the critical coupling at zero quark 
mass, 6/g 2 (0), we use the scaling relation for the 
crossover (pseudocritical) coupling 

6/3p C = 6/ 5 2 (0) + x^Ntamg) 1 ^ 6 

which works well. 

Neglecting nonscaling contributions to ipip, we 
construct the QCD critical scaling function (for 
x > 0) using 0(4) critical exponents || and com- 
pare with the scaling function for 0(4) fijl in the 
figure below. Vertical and horizontal displace- 
ments of the log-log scaling curves are permitted 
and correspond to adjusting c x and c y . 
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The line segments locate the crossover in the 
QCD data and the dashed line on the 0(4) curve, 
the crossover in 0(4). The newer data are plotted 
with octagons and squares. 

We see: (1) The new data show generally worse 
agreement with the 0(4) scaling curve. (2) In 
terms of these scaling variables and at these quark 
masses, the 0(4) critical region is quite small, if it 
is to be found at all. Equivalently, nonnegligible 
nonscaling contributions may be masking critical 
behavior in this regime. (3) Nonscaling contribu- 
tions could be displacing the crossover location. 

4. LARGER N t 

A similar analysis at larger Nt is shown below. 
Perhaps there is improvement with increasing Nt. 
However, our N t = 12 sample includes data only 
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at a single quark mass, making it our weakest 
example. Furthermore, for Nt > 4 we have no 

results for L > 2N t . 
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5. SPECULATIONS 

New results at smaller quark mass and larger 
volume at Nt = 4 have raised doubts about the 
extent of the previously observed agreement be- 
tween QCD and 0(4) In striking contrast 
with Wilson fcrmions with an improved gauge ac- 
tion Q , the conventional staggered fermion action 
with the conventional choice of scaling variables 
does not show good agreement with the 0(4) scal- 
ing function at present quark masses and temper- 
atures. 

Why? (1) Perhaps the critical region is at- 
tained only when tt and a correlation lengths are 
considerably greater than 1/T C . Here, typically, 
they are smaller. (2) At strong coupling the vari- 
ables 6/g 2 and ma may have a strongly nonlinear 
correspondence with spin model temperature and 
magnetic field T, h, which becomes approximately 
linear only in a small region around the critical 
point. The nonlinearities may become weaker in 
the weak coupling limit and perhaps may be di- 
minished with improved actions. Perhaps with a 
more suitable choice of scaling variables, we may 
find better agreement. 
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